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ABSTRACT 

We study the ultraviolet to far-infrared (hereafter UV-to-IR) SEDs of a sample 
of intermediate redshift (0.2 < z < 0.7) UV-selected galaxies from the ELAIS-Nl 
and ELAIS-N2 fields by fitting a multi- wavelength dataset to a library of GRASIL 
templates. Star formation related properties of the galaxies are derived from the 
library of models by using the Bayesian statistics. We find a decreasing presence 
of galaxies with low attenuation and low total luminosity as redshift decreases, 
which does not hold for high total luminosity galaxies. In addition the dust 
attenuation of low mass galaxies increases as redshift decreases, and this trend 
seems to disappear for galaxies with M* > lO^-*^ Mq. This result is consistent with 
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a mass dependent evolution of the dust to gas ratio, which could be driven by a 
mass dependent efficiency of star formation in star forming galaxies. The specific 
star formation rates (SSFR) decrease with increasing stellar mass at all redshifts, 
and for a given stellar mass the SSFR decreases with decreasing redshift. The 
differences in the slope of the M*-SSFR relation found between this work and 
others at similar redshift could be explained by the adopted selection criteria of 
the samples which, for a UV selected sample, favours blue, star forming galaxies. 

Subject headings: surveys: GALEX, SWIRE — ultraviolet: galaxies — galaxies: 
evolution 



1. Introduction 

Multi-band observations of galaxies are becoming increasingly important for a coherent 
understanding of the formation and evolution of galaxies. The increasing availability of wide 
field surveys at different wavelengths (e.g. GALEX, Spitzer, 2MASS, SDSS, 2dF) provides 
a wealth of data to better constrain the evolution of galaxies by analyzing it from several 
physical standpoints. In the last decade, wide area surveys have evidenced a strong evolution 
of some galaxy properties with redshift, particularly those related to star formation and dust 
attenuation. In particular, the local star formation rate density was found to be far lower 
than it was at z ^ 1 (Lilly et al. 1996; Madau et al. 1996; Flores ct al. 1999; Hopkins 
2004), although there is no clear agreement on the epoch where it was maximum (Heavens 
et al. 2004). A strong evolution of luminosities with redshift is also suggested by the study 
of the rest frame ultraviolet (UV; Arnouts et al. 2005) and rest frame infrared (IR; Xu 2000; 
Chary & Elbaz 2001; Le Floc'h et al. 2005) luminosity functions at different redshifts. 

The overall dust attenuation in star forming galaxies also shows signatures of evolution 
with redshift. Takeuchi et al. (2005a) have reported a strong decrease in the ratio of the 
IR to UV star formation rate (SFR) densities {pir/ puv) between z — 1 and z — 0, which 
can be interpreted as a decrease over this redshift range in the dust attenuation of galaxies 
averaged over a volume element. The situation is far from clear when individual galaxies are 
studied. Buat et al. (2007) and Xu et al. (2007) report a slightly lower dust attenuation 
for LlRGs at 2; ~ 0.6 — 0.7 compared to LIRGs aX z = which seems not to be true for 
ULIRGs (Xu et al. 2007). Burgarella et al. (2007) reported a lower dust attenuation (as 
measured from the IR/UV ratio) for LBGs at 2; ~ 1, compared to the sample of UV selected 
local galaxies of Buat et al. (2006), and these LBGs are more attenuated than the sample 
of Reddy et al. (2006) at ^ ~ 2. 
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Another quantity which has been extensively studied is the specific SFR (SSFR), which 
is defined as the ratio of the SFR to the stellar mass of the galaxy. This SSFR is now used 
to study the downsizing phenomenon - i.e. when z decreases intense star formation activity 
becomes limited to less and less massive galaxies - from z ~ 5 to present time (e.g. Cowie 
et al. 1996; Heavens et al. 2004; Brinchmann & Ellis 2000; Boselli et al. 2001; Bauer et al. 
2005; Feulner et al. 2005). High mass galaxies with values of SSFR as high as those of low 
mass galaxies have been reported in the local Universe, most of them being LIRGs and/or 
ULIRGs (Hammer et al. 2005; Iglesias-Paramo et al. 2006; Buat et al. 2006), but they are 
far less numerous than the massive galaxies with negligible star formation activity. Another 
interesting result reported by Bell et al. (2005) is that the SSFR of high mass galaxies was 
much higher at z ~ 0.7 than in the local Universe and could be responsible for the steep 
decrease in the SFR density from this epoch to the present one. 

Synthetic templates of star forming galaxies have demonstrated to be powerful tools to 
derive star formation related properties of galaxies, especially in the wavelength range UV-to- 
NIR (e.g. Leitherer & Heckman 1999; Bruzual & Chariot 2003; Fioc & Rocca-Volmerange 
1997). Furthermore, many efforts have been devoted to model the properties of the dust 
emission in star forming galaxies (e.g. Dale et al. 2001; Chary & Elbaz 2001; Rowan- 
Robinson et al. 2004). There have also been several attempts to reproduce the whole 
SEDs of star forming galaxies by using different methods: Silva et al. (1998) and Piovan 
et al. (2006) take into account the effects of the interstellar medium by considering a two 
phase model with a diffuse component and a dense one corresponding to molecular clouds. 
Devriendt et al. (1999) synthesize the stellar component and use a phenomenological fit 
to compute the extinction curve assuming a simple geometry for the dust, which yields the 
fraction of the stellar light reprocessed by the dust. Popescu et al. (2000) solve the radiation 
transfer problem for a finite disk and a bulge by using a consistent treatment of grain heating 
and emission. Finally, Dopita ct al. (2005) perform a detailed modeling of the PAH emission 
features and of the reprocessing by the dust grains over a range of sizes. 

In this paper we derive properties of a sample of UV selected galaxies by fitting multi- 
band data to a set of GRASIL (Silva et al. 1998) templates in a statistical way, similar 
to what has been done for other samples of galaxies by using UV-to-optical templates (e.g. 
Kauffmann et al. 2003) or UV-to-optical templates plus a total IR luminosity (Burgarella 
et al. 2005) built from other population synthesis codes. We obtain mean star formation 
related properties at different redshift bins and explore the likely evolution of these properties 
with redshift. The paper is organized as follows: In § 2 the multi-band dataset is described. 
The main properties of the synthetic models are explained in § 3. The fitting method and 
the determination of the parameters are described in § 4 and 5 respectively. § 6 contains a 
discussion on the main results obtained, and a brief summary is presented in § 7. Throughout 
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this paper we will use the following cosmological parameters: Hq = 72 km s ^ Mpc ^, 
Qm = 0.3 and Qx = 0.7. 

2. The dataset 

2.1. Sample selection 

Our dataset contains multi-wavelength (UV to FIR) photometric data of galaxies in the 
ELAIS fields. The UV data were taken by GALEX (Martin et al. 2005a) and consist of 19 
circular frames (0.6 deg radius) in the ELAISNl and ELAISN2 fields, observed through two 
filters with the bandpass centered at 1530A (FUV) and 2315A (NUV). In this paper we use 
only NUV photometry, since FUV photometry is available only for a few of these frames. 
The final frames were the result of the co-add of different observations with typical total 
exposure times of ~ 15000 s. The NUV photometry used in this paper corresponds to the 
GALEX Second Data Release (available at l| http : / /galex. stsci . edu/ G R2 / [ ) . Table [1] shows the 
coordinates and exposure times of the GALEX frames used in this work. 

The optical to FIR data correspond to the Spitzer Wide-area InfraRed Extragalac- 
tic (SWIRE) survey and consist of 12-band photometry: U,g',r',i', Z taken at INT2.5m 
with the WFC; 3.6/im,4.5/im,5.8/im,8.0/im from Spitzer/IRAC; and 24/xm,70/im,160/im from 
Spitzer/MIPS. In addition, photometric redshifts are provided for all the galaxies on the list. 
Throughout this paper, the merged catalog (version June 2005, Rowan-Robinson, private 
communication) was used. Details on the photometry of the SWIRE sources and on the 
quality of the photometric redshifts are given in Rowan- Robinson et al. (2005, hereafter 
RR05). 

This work intends to study the properties of a UV selected sample of galaxies so we 
focused on the GALEX frames and selected all the sources with mjvc/y < 22.5 mag. Sources 
located out of the inner 0.5 degree radius from the center of the GALEX frames were rejected 
since going to larger distances from the center can introduce spurious detections and ghosts 
that often appear close to the border. The angular resolution of the GALEX frames (FWHM 
~ 5") does not allow a proper star/galaxy separation, thus a cross-correlation between the 
UV and SWIRE catalogs was performed. All the GALEX sources for which a single SWIRE 
source (classified galaxy in the g' and r' bands) is present in a circle of a 2" radius were 
included in our sample. In order to avoid confused detections due to the different angular 
resolution of the GALEX and (SWIRE) optical frames, those UV sources for which more than 
one optical detection was present within a circle of radius 6" were excluded from the final 
sample. In this way, a magnitude limited UV selected sample of galaxies free of confusion 
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effects was built. 

Both, the 6" isolation condition previously mentioned and the cross-correlation with the 
SWIRE catalog may prevent a pure UV selection of our sample due to two reasons: (a) We 
may lose closely interacting galaxies and mergers which are either not resolved in the NUV 
frames (but they do in the SWIRE catalog) or are so close that they break the 6" isolation 
condition; and (b) The SWIRE catalog contains only sources which were detected at both 
r' and 3.6/im, which means that our final sample could lack blue galaxies included in our 
original NUV list but not detected either at r' or at 3.6/im. However, the completeness limits 
at 3.6/im (3.7/iJy) and r' (23 Vega mag) are much deeper than our NUV survey so including 
only SWIRE sources classified as galaxies does not introduce any bias in the sample (see 
Xu et al. 2007 for a discussion on this point). Also shown in Figured] are the NUV counts 
derived from our sample compared to those of Xu et al. (2005). As can be seen, both sets 
of counts are consistent within the errors. Thus, we conclude that the selection procedure 
does not include significant biases and for practical reasons, our UV selected sample can be 
considered complete down to rriNuv = 22.5 ABmag. 

In the end we selected 6980 galaxies with (photometric) redshifts in the interval 0.2 < 
z < 0.7. There are several reasons for this redshift restriction: firstly, the relative uncertainty 
of the photometric redshifts increases as they approach (see RR05 for details), which will 
surely induce large uncertainties in the derived parameters dependent on redshift; secondly, 
Buat et al. (2006) have already analyzed the properties of a UV selected sample of galaxies at 
z «i 0; finally, going to higher values of redshift would increase the number of non detections, 
mainly at IR wavelengths, and this fact would hamper our analysis. Figure [2] shows the 
redshift distribution of our sample galaxies and Table [2] the fractions of detected galaxies at 
all the available bands. 



2.2. AGN contamination 

Although a UV selected sample of galaxies is expected to be dominated by star forming 
galaxies, we cannot neglect the possibility of some contamination due to the presence of 
AGNs. In order to estimate the fraction of AGNs, we use the classification of RR05. In 
that papet the authors fitted the data usmg a set ot 8 empirical optical-MIR template^ 
corresponding to E, Sab, Sbc, Scd, Sdm, Sburst, AGNl and AGN2. Figure [3] shows the 
distribution of their best fit templates for our 6980 galaxies. As can be seen, only about 
0.4% of the galaxies were classified as AGN by RR05. In addition, less than 4% of the 



^No UV fluxes were used in RR05. 
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galaxies were classified as E, the rest of tliem being classified as spirals or starbursts. This 
result supports the fact that from the optical point of view our sample is dominated by star 
forming galaxies, with very low contamination from AGNs and ellipticals. 

3. Synthetic SEDs 

A set of synthetic templates of SEDs were computed with the code GRASIlB. This code 
treats with particular care the effects of a dusty interstellar medium for different geometries 
of the dust and stars and is coupled with a chemical evolution code (CHE_EVO) that takes 
into account in a consistent way the chemical evolution of the galaxy and the star formation 
history. This code describes one-zone open models depending only on time, including the 
infall of primordial gas. A more detailed description of this code can be found in Silva et al. 
(1998). Although GRASIL allows the variation of a large set of free parameters (see Silva et 
al. 1998 for a detailed description of the code), to build our grid of templates only a small 
subset of parameters were varied (Silva, private communication). The values spanned by the 
free parameters were chosen so that they cover the values proposed by Silva et al. (1998) to 
reproduce the SEDs of star forming galaxies of different Hubble types in the local Universe 
(see Silva et al. 1998 for details on fits to spiral and starburst galaxies). 

• The SFR was assumed to follow a Schmidt-type law of the form 

SFR{t) = vMgitf (1) 

where Mg{t) is the gas mass at any time, k = 1, and u is a free parameter varying from 
^ 2.56 X 10"^ to ^ 4. 

• Gas is continuously infalling at a rate proportional to exp(— t/rj„/), where Tinf is a free 
parameter taking values from 0.1 to 21.6 Gyr. These values extend above and below 
the typical values found for the Milky Way by Boissier & Prantzos (1999): r = 1 Gyr 
at i? = 1 kpc, and r = 10 Gyr at i? = 18 kpc. GRASIL does not allow us to vary 
infall rates with the distance to the center of the galaxy, so single values are adopted 
for each model. 

• GRASIL also allows the addition at any time of an analytical extra star formation law 
(simulating a burst of star formation) combined with a Schmidt law. We chose an 



^Details about the use of the code can be found at the web site 
http://web.pd.astro.it/granato/grasil/grasil. html , A web interface to use GRASIL can be found at 
http://web.pd. astro. it/galsynth/index.php[ 
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exponentially declining functional form with declining timescale equal to 50 Myr. For 
comparison, Kauffmann et al. (2003) and Salim et al. (2005) use bursts of constant 
intensity between 30 and 300 Myr after the beginning of the burst. Four intensities of 
the burst were allowed in our work, corresponding to 0.33%, 1%, 3% and 9% of the 
total mass infall. 

• GRASIL assumes two components of the interstellar medium (ISM): the molecular 
clouds - containing the young stars and a dense ISM -, and the diffuse gas (cirrus) 
- containing stars that escaped from the molecular clouds (free stars) and a diffuse 
ISM -. The typical time required to escape the molecular clouds, tesc, is also a free 
parameter that we let vary from 1 to 4 x 10'' yr, a value consistent with the one used 
by Chariot & Fall (2000) and Panuzzo et al. (2007). An escape time of 2 x 10^ yr is 
able to reproduce the SEDs of some local spiral galaxies (Silva et al. 1998). 

• The optical depth of the molecular clouds, tmc oc mMc/f^MCj where mMc and vmc 
are the mass and radii of the molecular clouds respectively, was taken as another free 
parameter in our set of templates. We fixed vmc = 14 pc and let Mmc vary from 
3 X 10^ to 3 X 10® Mq. These values are close to the ones measured by Rosolowsky et 
al. (2007) for Giant Molecular Clouds in M33, although their masses are slightly lower 
and their radii slightly larger than the ones proposed here. On the other hand, our 
values are close to the ones used in Silva et al. (1998) to reproduce a sample of local 
spirals. 

• GRASIL allows us to stop the evolution of the galaxy at any time, so that galaxies of 
different ages can be simulated. We assumed that all galaxies are coeval with an age 
of ~ 13 Gyr. For the redshift interval selected in this work, and given that SEDs are 
almost insensitive to small age differences, we produced models of 8 and 12 Gyr. 

• The inclination angle of the galaxies may play a role in the observed SED of a galaxy, 
since the more inclined a galaxy is, the more attenuated it appears at UV and optical 
wavelength^. This may result in an incorrect derivation of the dust attenuation and 
star formation related parameters. Although we lack available information on the 
inclination angles for our galaxies, we generated models corresponding to inclination 
angles of 15°, 45° and 75°. 

In all our models, stars and gas are distributed following exponential disks with scale 
length 4 kpc and scale height 0.4 kpc. These values are of the order of the ones reported for 



•^The IR emission remains almost unchanged since it is isotropic 
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large samples of spiral galaxies (de Jong et al. 1996; Yoachim & Dalcanton 2006) although 
other authors report slightly higher values for the scale height (Barteldrees & Dettmar 1994). 
The assumed fraction of molecular gas is 50%, which is rather high but still reasonable for 
star forming galaxies (Boselli et al. 2002). Finally, the dust to gas ratio is assumed to be 
proportional to the metallicity and normalized to the Milky Way value at solar metallicity. 

Our database is composed of a total of 5127 templates. Table [3] lists the free parameters 
of the templates used in this paper. 

The evaluation of the total IR luminosity between 8 and lOOOyum {Lju) is a potential 
source of uncertainty, mainly when almost no detections at FIR wavelengths and detections 
at 24yum for only about half of the sources are available, as is the case for our dataset. For 
this reason, we prefer to be consistent with previous observational results and we restrict 
our models to those which satisfy the following relations: 

Ljn = 10°-2« X L^-o^° (2) 

Ljn = W'-'' X L°f « (3) 

L,n = 10°-^° X Leo (4) 

where Lg, L25 and Lgo are the monochromatic luminosities at 8, 25 and 60/im respec- 
tively. These three observational constraints have been reported by Caputi et al. (2007) 
and Takeuchi et al. (2005b) and hold over several orders of magnitude. In order to account 
for the dispersions found around these relations, we allow models to deviate by a factor of 
3 above and below these average relations. We note that whereas two of these relations are 
found for local galaxies (Takeuchi et al. 2005b), the one by Caputi et al. (2007) holds for 
galaxies at < 2; < 2. Thus, by applying these local relations to our models we assume that 
the IR SEDs are suitable to estimate the total IR luminosities of higher redshift galaxies, as 
seems to be the case (Egami et al. 2004; Le Floch et al. 2005; Marcillac et al. 2006). 

In order to show how our models reproduce the SEDs of the sample galaxies from the 
UV up to the IR, in Figures H] to [H] several observed color-color plots for the galaxies and 
for the templates at different redshifts are presented. The observed colors involving only 
optical bands (figure H]) are properly reproduced by the models at any redshift. The fraction 
of galaxies located at distances larger than 0.2 mag from any template is lower than 3% at 
all redshifts. Concerning the observed colors involving optical and IR data (figure [5]) the 
situation is quite similar, with less than 4% of the galaxies more distant than 0.2 mag to 
any template. However, the situation is less favourable for the observed colors involving 
the NUV fluxes, as is shown in figure [61 In this case, the fraction of galaxies showing 
{NUV — g') colors not reproduced by the models to an accuracy of 0.2 mag increases from 
5% for z G [0.2,0.3] to 25% for z G [0.6,0.7]. The reason for this discrepancy between data 
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and models is not clear. A careful revision of the most extreme cases yields NUV fluxes 
much larger (compared to their optical fluxes) than those expected for star forming galaxies. 
The median NUV fluxes of the discrepant galaxies are similar to that of the total sample, 
so photometric uncertainties can be discarded as the only reason for this effect. In any case, 
the total number of discrepant galaxies is less than 10%, therefore most of the galaxies can 
be properly fitted by the models. 



4. SED fitting procedure 

For each galaxy of our sample the reduced corresponding to all the available templates 
satisfying Eqs. [21 [3] and H] is computed. In this way the scaling factors for all the templates 
are obtained. Many of our galaxies have no available IR fluxes, so the upper limits in the 
derivation of are included. 

Following the usual notation, is defined as: 

i 

where n is the number of bands at which the galaxy was detected or for which the model 
flux exceeds the upper limit flux at any wavelength, and 

{{f\,i - JKifl {^f\A? , if detection at \i 
- Jx,i)V{Aflf , if non detection at A, and fl, < jx,i 
, if non detection at Aj and f^i > 

(6) 

where f\^i is the detected flux at Aj, A/a^ is the flux uncertainty at Aj, f^^ is the upper limit 
at Aj, Af^^ is the uncertainty associated to f^i, and jA,i is the model flux at Aj. 

For each galaxy, the associated value of the reduced is the minimum among all 
the templates. Figure [7] shows the histogram of the reduced the median value being 
2.58. In order to illustrate the quality of the fits for different values of a selection 
of them are shown in Figure [H It can be seen that the quality of fits with x^ < 4 is 
reasonably good. For our analysis, we prefer to keep only galaxies with reasonably low 
values of x^ since otherwise, the uncertainties in the derived parameters will be large. For 
this reason, we exclude all the galaxies for which x^ > 4, and we end up with a sample of 
5133 galaxies. We have verified that discarding these galaxies does not introduce any bias in 
our sample: the fraction of galaxies detected at 24/im, the redshift distribution of galaxies 
and the average values of rest frame Lpjjy are similar for the complete and the restricted 



(good quality) samples. In addition, by imposing a limiting we also reduce the fraction 
of galaxies with discrepant {NUV — g') color which is finally of 1%, 1.5%, 7%, 5% and 11% 
for z e [0.2,0.3], [0.3,0.4], [0.4,0.5], [0.5,0.6] and [0.6,0.7] respectively. 



5. Determination of parameters and associated uncertainties 

For each template, GRASIL provides a complete SED (UV to FIR), a list of parameters 
related to the evolution of the galaxy such as the SFR and mass assembly as a function of 
time. For our purposes, we derive the luminositjill at any wavelength from the output SED. 
The total Ljpi is the luminosity emitted between 8 and 1000/im. The SFR and stellar mass 
(M*) are directly given by GRASIL. In this paper we will consider the SFR averaged over 
the last 10^ yr, since this is the typical timescale of star formation to which the UV and IR 
fluxes are sensitive (e.g. Iglesias-Paramo et al. 2006). 

The physical parameters associated to a given galaxy are the median values of the 
likelihood distribution resulting from applying the bayesian statistics to the whole set of 
templates, as described in detail in Kauffmann et al. (2003). The confidence intervals are 
defined as those containing 95% of the likelihood distribution. 

Figures P to [TT] show the distribution of the confidence intervals corresponding to Lpuv-, 
M* and Ljr- As can be seen, the Lpuv and M* distributions are clearly unimodal. This 
is not unexpected since the contribution of the stars to the SED mostly ranges from the 
UV to the NIR, which is reasonably well covered in our dataset. The median values of the 
confidence intervals are 0.34 and 0.52 dex respectively for Lpuy and M*. 

On the contrary, the Lm distribution shows a bimodal shape, with a median value of the 
confidence interval of 0.73 dex. This different behavior is due to the paucity of the coverage 
at FIR wavelengths of our dataset. In fact, the fraction of galaxies detected at 24/70/160yum 
(the wavelengths contained in the interval along which Lir is computed) is about 47/3/l9fl. 
This means that for about 50% of the sample, Lm is extrapolated from the GRASIL SEDs 
which are constrained mostly at UV-optical-MIR wavelengths. 

In order to evaluate the quality of the values of Lir derived from the GRASIL models, 
we compared our results with the models of Chary & Elbaz (2001). These models provide 
whenever the redshift of the galaxy and the flux at 24/im are known. For this reason. 



^Throughout this paper, monochromatic luminosities are defined as \F\, where A is the reference wave- 
length and F\ is the flux density at this wavelength. 

^These fractions correspond to the good quality sample after excluding the galaxies with > 4. 
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a direct comparison is only possible for the subsample of galaxies detected at 24/xm. The 
comparison is shown in Figure [121 Our determination of Lm is larger than that of Chary 
& Elbaz (2001) by 25% on average, with a dispersion of 0.13 dex. If we restrict ourselves 
to those galaxies detected at 70 /im, our estimation is still 4% larger with a dispersion of 
0.12 dex. 



6. Results and discussion 

6.1. Dust attenuation versus total luminosity and stellar mass 

Lboi (defined as Lpuv + Lir) is a good tracer of the total luminosity emitted by young 
stars since it contains both the contribution of the light directly emitted by young stars plus 
the IR emission of the dust, which absorbs and scatters part of this light (e.g. Martin et al. 
2005b; Buat et al. 2006). The energetic budget between the IR and UV emission (namely 
the Lm/ Lpuv ratio) is known to be a good tracer of the dust attenuation at UV wavelengths 
(Xu & Buat 1995; Meurer et al. 1999; Gordon et al. 2000). It has already been reported 
in the literature that galaxies with large SFRs also tend to show a larger dust attenuation 
(Wang & Heckman 1996; Buat & Burgarella 1998; Hopkins et al. 2001). Similar relations 
have also been found for galaxies at higher redshifts (Burgarella et al. 2006; Reddy et al. 
2006) selected at UV rest frame. With our sample we have the opportunity to study the 
evolution with redshift of this relation for UV selected galaxies from z ^ 0.7 to z ~ 0.2. 

Figure fT3l shows Lpuv/Lm as a function of Lboi- We have split our sample into five 
redshift bins: [0.2,0.3], [0.3,0.4], [0.4,0.5], [0.5,0.6] and [0.6,0.7]. The top-left panel also 
shows the average relation of Buat et al. (2006) for UV selected local galaxies with the same 
limiting magnitude as our subsample with z G [0.2,0.3]. As can be seen, for Lhoi > 10^° 
the agreement is quite remarkable with the data at 2 G [0.2,0.3]. For the sake of simplicity 
and homogeneity we constraint our study of the evolution of the Ljji/Lpuv ratio to our 
subsamples at different redshift bins. 

The comparison between the different redshift bins is difficult because the observational 
limit [mjsfuv ^ 22.5 ABmag) translates to different luminosity limits, thus we can only 
samples of galaxies in regions of the plot which are free of observational biases can be 
compared. The observational limits for each redshift bin were estimated by assuming that 
the limiting NUV magnitude translates into Lpuv rest frame and obtained limiting values of 
Lpuv = 10^ 2.25 X 10^ 4.95 x 10^ 8.46 x 10^ and 1.32 x 10^° Lq for z = 0.2, 0.3, 0.4, 0.5 and 
0.6 respectively. The limit at ^ = 0.6 restricts the comparison to a very reduced number of 
galaxies at lower redshifts in Figure [131 For this reason we perform our comparison within 
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regions (a) and (b) of the figure, which are free of observational biases for galaxies with 
z G [0.2, 0.6]. These two regions correspond to L^/ G [2.5 x 10^°, 10^^] and [10i\ 4 x 10^^] L©. 

We rely on the fact that galaxies should be equally distributed in any region of the 
diagram free of observational biases, if the dust attenuation does not evolve with redshift for 
a given total luminosity. Thus, a detailed study of the distributions of galaxies in regions 
(a) and (b) is performed at different redshift bins. The results of a Kolmogorov-Smirnov 
test comparing the populations of both regions for different redshifts are shown in Table |H 
For region (a) only galaxies with z G [0.3, 0.4] show distributions consistent with those at 
z G [0.2, 0.3] at the 95% confidence level. For the rest of the redshift bins, the distributions 
of galaxies are significantly different to each other. On the other hand, the distributions 
of galaxies in region (b), are consistent with being drawn from the same population at the 
95% confidence level, so no evolution with redshift is reported for them. Figure fT^ a) shows 
that the histograms of Lm/Lpuv for region (a) are shifted towards higher values as redshift 
decreases, which suggests that low total luminosity galaxies - region (a) - with low values of 
LiFt/Lpuv become scarce at low values of z. This decreasing presence of galaxies with low 
Lin/Lpuv as redshift decreases is not seen for high total luminosity galaxies - region (b) - 
as the Kolmogorov-Smirnov test states and Figure [TWb) suggests. 

The evolution of the dust attenuation with the stellar mass is also studied. Figure [T5] 
shows M* vs. Ljr/Lfuv at different redshifts is also studied. Again we are only limited by 
our UV selection because the rest of the relevant parameters (M*, Ljpi) are derived from 
the templates after fitting to the observational data. The only limit in Figure [TS] at a given 
redshift is a minimum M* associated to the minimum Lpuv imposed by the UV selection. 
Thus, the distribution of stellar masses represented in this figure corresponds to that of a 
UV selected sample and does not intend to be complete down to any value of M*. 

In Figure dn] the volume weighted average M* vs. Ljji/Lpuy at the different redshifts 
is shown. The average and standard deviation (corresponding to the error bars in the figure) 
values have been estimated as in Buat et al. (2006): 

{^ir/^^fuv) ^; [I) 

and 

2/r /r ^ ^i'^iiiLm/ LFuv)i - {Lm/ Lpuv)) 

[Lir/ Lpuv) — [o) 

where Ui = l/Vmax for galaxy i. 

It can be seen that Lj^/ Lpuy increases with stellar mass at all redshifts for galaxies with 
M* < 10^^ M0. A galaxy with M* ^ 2 x 10^° M© shows Ljr/Lfuv a factor of ~ 3 higher 
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at 2 = 0.65 than at z = 0.25. This difference is lower for a galaxy with M* ^ 7 x 10^° M0, 
a factor of ~ 2. For galaxies with M* > 10^^ Mq the trend seems to flatten and Lj^/Lpuv 
reaches an almost constant value with redshift. The plot also shows that for galaxies more 
massive than logM* = 11.5 the evolutionary trend of Li^/ Lpuy seems to invert in the sense 
that it decreases with redshift. Whether this effect is real or not is still to be confirmed since 
the statistics are scarce for this last mass bin. 

A tentative explanation of this evolution of Ljji/Lpuv with the stellar mass could be 
associated with the total dust content of star forming galaxies, which is linked to the metal- 
licity. A mass-metallicity relation has been reported at different redshifts in the sense that 
more massive galaxies are also more metal rich (Zaritsky et al. 1994, Tremonti et al. 2004). 
We also find that at a fixed stellar mass, the dust attenuation increases as redshift decreases. 
This result is consistent with the findings of Savaglio et al. (2005) and Erb et al. (2006) 
that for a given stellar mass, galaxies were less metallic at high redshift than they are at 
present time. This result suggests a mass dependent rate of chemical enrichment for galaxies 
that can be explained by invoking a mass dependent star formation history. In fact, N-body 
simulations by Brooks et al. (2007) claim that low star formation efficiencies are required 
for low mass galaxies in order to explain the mass-metallicity relation. We can estimate the 
star formation efficiency of our galaxies through the fitted parameters of GRASIL. The u 
parameter is the main driver of the SFR when we assume a Schmidt-type law (Eq. 1). Low 
values of u allow a fairly fiat SFH until the present epoch. On the contrary, high values of u 
imply a very intense star formation activity in the first epochs of a galaxy followed by a steep 
decrease of the SFR due to a quick exhaustion of the gas reservoir. This decrease in the SFR 
of massive galaxies has already been reported by Bell et al. (2005), Juneau et al. (2006) and 
Pannella et al. (2007) for different samples of galaxies from z ^ 2 to present time. Figure [T71 
shows the relation between M* and u for our galaxies. It is apparent that low mass galaxies 
show lower values of u while the opposite holds for high mass galaxies. With this downsizing 
picture in mind we can explain the differential evolution of dust attenuation as partially 
linked to the SFH of galaxies: high mass galaxies experience a first epoch of intense star 
formation activity, thus they will efficiently increase their metal (and dust) content. As the 
SFR begins to decline up to the present epoch the metal enrichment becomes slower as does 
the dust production. Conversely, low mass galaxies show almost constant star formation 
efficiencies, which translate into moderate but continuously increasing chemical enrichment 
and dust production. To illustrate this hypothesis the volume weighted metallicity evolution 
for different bins of M* is shown in Figure [181 Metallicities have been computed by the 
chemical evolutionary code CHE_EVO coupled to GRASIL in a consistent way with the rest 
of the parameters, and we remind the reader that in our models the dust to gas fraction 
was assumed to be proportional to metallicity. As can be seen, whereas for galaxies with 
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logM* G [11.5, 12] the metallicity increases by a factor of ~ 1.6 from z = 0.65 to z = 0.25, 
for galaxies with logM* G [10,10.5] the increasing factor is ~ 8. Thus the evolution of 
Liji/Lpjjy observed in figure [TBI could be understood in terms of an evolution of metallicity. 

6.2. The specific SFR 

The specific star formation rate (SSFR) is defined as SSFR = SFR/M*, that is, the 
ratio of the current SFR to the total stellar mass produced throughout the galaxy's lifetime. 
It is related to the birthrate parameter b (defined as SFR/ (SFR), where (SFR) is the 
averaged SFR throughout the galaxy's lifetime), but it does not depend on the timescale 
required to build all the stellar mass. If we assume that the SFR of the galaxy will remain 
almost constant for a long period of time, the SSFR is equivalent to the inverse of the time 
required by a galaxy to double its stellar mass. 

Figure [T^ shows the SSFR as a function of the stellar mass. The observational limits are 
indicated with dashed lines. The limits correspond to the minimum stellar mass and SFR 
allowed for each Lpuv of our sample. The horizontal lines in each panel indicate the SSFR 
required for a galaxy to double its stellar mass from each epoch to present time assuming 
constant SFR. It is found that less massive galaxies tend to have a larger SSFR, consistent 
with the scenario already reported by several authors (Cowie et al. 1996; Brinchmann & 
EUis 2000; BoseUi et al. 2001; Brinchmann et al. 2004; Heavens et al. 2004; Feulner et al. 
2005; Bell et al. 2005; Buat et al. 2006) for samples of galaxies selected on the basis of 
different wavelengths and at different redshifts. For z G [0.2, 0.3] most of the galaxies are 
below the doubling-mass line, but this trend changes at higher redshifts and at 2; G [0.6, 0.7] 
most of the galaxies are well above thi= linfi TM. re.ult is clearer in Fi gure [20] where the 
volume weighted average SSFR as a function of M* for the different redshift bins are shown. 

For a given M*, the SSFR is found to decrease with redshift. The slope of the variation 
of the SSFR with M* as well as the rate of decline of SSFR with z for different stellar masses 
are clues to understand the origin of the downsizing process. We find a shallower variation of 
SSFR versus M* (i.e. a larger SSFR for massive galaxies) as compared to other works aimed 
at studying the evolution of the SSFR-M* variation and covering (at least) the same redshift 
range (e.g. Feulner et al. 2005, Bauer et al. 2005, Zheng et al. 2007). From Figure [20l the 
SSFR is also found to decrease by a larger factor for massive galaxies than for low mass ones 
from z = 0.7 to z = 0. As an example, the SSFR lowers by a factor of ~ 8 from z ~ 0.65 



^Feulner et al. (2005) call the galaxies falling below and above this line "quiescent" and "star-forming' 
respectively. 
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to 0.25 for galaxies more massive than 10^^ M©, whereas this factor is of the order of 4 for 
galaxies with M* 10^'^ Mq. This result differs from the recent findings of Zheng et al. 
(2007) who conclude a similar decrease of SSFR with redshift for all stellar masses from 
z = 0.9 to present time. 

In this respect, the wavelength at which samples are selected may introduce some bias 
on the obtained results. To illustrate this point, Figure [21] shows M* vs. rest frame {U — V) 
for our galaxies. The lines separating the red and blue sequences at the different redshifts 
(from Bell et al. 2007) are shown in separated panels. As can be seen, the majority of 
our galaxies belong to the blue sequence at all redshifts and almost no signs of any color 
bimodality are found (e.g. Strateva et al. 2001). This is a consequence of the UV selection 
applied to our sample. Contrary to that, the sample of Elbaz et al. (2007) selected at rest 
frame optical wavelengths shows a clear bimodality in the color magnitude plot. The same 
is found by Bell et al. (2007) for their sample selected at optical magnitudes. These authors 
also show that red galaxies show a steeper relation between M* and SSFR than blue ones. 
From a general point of view, care should be taken when comparing star formation related 
properties from samples selected with different criteria since the results will depend on the 
fraction of red to blue galaxies of each sample. In particular, our results show that UV 
selected samples are extremely efficient at selecting blue galaxies. However, as has been 
shown in previous papers (Buat et al. 2005; Iglesias-Paramo et al. 2006), a complete view 
of the evolution of star forming galaxies would require a sample of galaxies containing both 
the bright (UV) and dark (IR) sides of star formation. 



7. Conclusions 

In this paper a study has been performed on the properties of a sample of UV selected 
galaxies at intermediate redshifts (0.2 < z < 0.7). The properties were inferred from the fit 
to a grid of synthetic models constructed with the code GRASIL through a Bayesian analysis. 
Most of the galaxies fit reasonably well to the models, which indicates that a Schmidt-type 
law combined with recent bursts of star formation properly reproduces the observed SEDs 
of galaxies at this redshift range. 

It was found that, whereas low luminosity galaxies with low attenuation were signifi- 
cantly more numerous at high redshifts, this result is not seen for high luminosity galaxies. 
In the same way, dust attenuation is found to increase with decreasing redshift for galaxies 
with low M*, but this trend seems to disappear for galaxies with logM* > 11. This result is 
consistent with a mass dependent evolution of the dust to gas ratio, and thus of metallicity, 
which could be driven by a mass dependent star formation history of galaxies: the most 
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massive galaxies have formed most of their stars by 2; ~ 0.7 and thus their metaUicity would 
not significantly change from z — 0.7 to ^ = 0.2, and likewise their dust content. The oppo- 
site is found to hold for low mass galaxies, which experience a larger increase in metaUicity 
over the same period, producing the reported evolution in the dust attenuation. 

The SSFR were found to decrease with increasing stellar mass, as has been known for 
many years (e.g. Cowie et al. 1996), and the trend holds at all redshifts explored in this 
work. Assuming constant SFRs, the fraction of galaxies able to double their stellar masses 
from their epoch to the present time decreases with rcdshift. At 2; > 0.6 most of the galaxies 
with Af* > 2 X 10^*^ M arc above the doubling-mass line whereas the opposite is true for 
such massive galaxies at 2; ~ 0.2. A shallower decrease of the SSFR with increasing stellar 
mass is found than reported in previously published studies in the same redshift range and 
based on visible-NIR selections. This difference is likely to be due to the sample selections: 
whereas samples selected at optical/IR wavelengths include galaxies belonging to both the 
blue and red sequences, our UV selection is dominated by blue (star forming) galaxies. 
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Tablc 1: Main properties of the GALEX frames: (1) Id. of the GALEX frames; (2) Right 
Ascension of the center of the frame, in J2000.0 equinox; (3) DecUnation of the center of the 
frame, in J2000.0 equinox; (4) Total exposure time of the co-added NUV frame. 
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Table 2: Detection rates at all bands for the sample galaxies. 
Band NUV U ? T Z 

Det. rate (%) 100.0 90.0 100.0 100.0 100.0 98.1 



Band 3.6/Lim 4.5/im 5.8/im 8.0/im 24/im 70/im 160/im 



Det. rate (%) 100.0 91.9 40.1 62.7 49.9 4.6 1.6 
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Table 3: Free parameters of the GRASIL models. 



SFR SL, SL+0.33%, SL+1%, SL+3%, SL+9% 

t,sc (Gyr) 0.01, 0.02, 0.04 

Mrnol (Mq) 3 X 10-\ 10^ 3 X 106 

r (Gyr) 0.1, 0.6, 3.6, 21.6 

logi/ -3.6, -2.9, -2.2, -1.5, -0.8, -0.1, 0.6 

Age (Gyr) 8, 12 

Incl. Angle 15°, 45°, 75° 



Table 4: Results of the Kolmogorov-Smirnov test for the distributions of Lir/Lfuv a-t dif- 
ferent redshifts. 
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16 18 20 22 24 

NUV (ABmag) 

Fig. 1. — NUV counts derived from our UV selected sample after matching with the optical- 
MIR SWIRE catalog (red). NUV counts derived by Xu et al. (2005) are also shown (blue). 
The histogram with the number of galaxies detected for each magnitude bin is also shown. 
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Fig. 2. — Redshift distribution of the sample galaxies. The numbers on top of the histogram 
correspond to the number of galaxies in each redshift bin. 
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Fig. 3. — Histogram of the best fit optical templates found by RR05 for our sample galaxies. 
The fractions of galaxies corresponding to each bin are indicated. 
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Fig. 4. — Observed {g'—r') vs. {r'—i') (ABmag) for the galaxies (red dots) and the templates 
(black dots). 
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Fig. 5. — Observed {g' — r') vs. (r' — /s.e) (ABmag) for the galaxies (red dots) and the 
templates (black dots). 
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Fig. 6. — Observed {NUV — g') vs. {g' — r') (ABmag) for the galaxies (red dots) and the 
templates (black dots). 
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Fig. 7. — Reduced -y^ distribution of the sample galaxies. The vertical line indicates the 
median value = 2.58. 



- 31 - 



Chi2 = 0.29 z = 0.28 




10' 10" 10^ 10^ 10' 

lambda (Angstroems) 



Fig. 8. — Best fits for several galaxies of tlie sample. Blue symbols correspond to detections. 
Red arrows correspond to upper limits. Error bars are indicated. 
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Fig. 8. — Continued. 
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Fig. 8. — Continued. 
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Fig. 9. — Histogram of the width of the 95% confidence interval of Lpuv for our sample 
galaxies. 
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Fig. 10. — Histogram of the width of the 95% confidence interval of M* for our sample 
galaxies. 
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Fig. 11. — Histogram of the width of the 95% confidence interval of Ljr for our sample 
galaxies. 
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Fig. 12. — LjR from Chary & Elbaz (2001) models compared to Lm derived from the 
GRASIL models. Blue (Red) points correspond to galaxies detected at 24/im (70/im). 
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Fig. 13. — Lhoi vs. Lir/Lfuv for the sample galaxies. Grey lines indicate the observational 
limits above which no galaxies of a given redshift bin are allowed. For each panel the 
corresponding observational limit is shown in bold type. The limiting FUV luminosities 
were estimated assuming that the limiting value of the observed mj^uy = 22.5 ABmag 
applies to the rest frame FUV luminosity at all redshifts. These limiting values are Lpuv = 
10^2.52 X 10^4.95 x 10^8.46 x 10^ and 1.32 x 10^° L© for z = 0.2,0.3,0.4,0.5 and 0.6 
respectively. The relation found for UV selected local galaxies by Buat et al. (2006) is 
shown in red in the top-left panel. The filled grey regions are the areas (free of observational 
biases) over which we have estimated the average Lm/Lpuv for the different redshift bins. 
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Fig. 14. — Distribution of Lm/Lpuy for regions (a) (upper plot) and (b) (lower plot) of 
Figure [13 
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Fig. 15. — M* vs. Lm/Lpuv for the sample galaxies. 
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Fig. 16. — Volume weighted average Li^/Lpuy vs. M* for the different redshift intervals. 
Error bars correspond to (volume weighted) standard deviations. The mass intervals corre- 
spond to: M* e [109■^ 10^°], [10^°, [10l°•^ 10"], [10", 10"-^], [10"■^ W\ The points 
have been slightly shifted horizontally in order to avoid confusion. 



-42 - 



Fig. 17. — M* vs. V for the sample galaxies. For each galaxy, the value of v is the average 
over all the models weighted by exp(— 
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Fig. 18. — Redshift vs. weighted averaged metalhcity (in solar units) for different stellar 
mass. 
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Fig. 19. — M* vs. SFRg/M* for the sample galaxies. The grey lines correspond to the 
(approximate) observational limits for each redshift bin. In each panel the observational 
limit corresponding to the indicated redshift bin is shown in bold type. The average relation 
found for UV selected local galaxies by Buat et al. (2006) is represented by a dashed cyan 
hne. The horizontal hnes represent the required SSFR to double a galaxy's mass between 
each redshift epoch and the present. 
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Fig. 20. — Volume weighted average values of SSFR as a function of M*. The mass intervals 
are: logM* G [9, 9.5], [9.5, 10], [10, 10.5], [10.5, 11], [11, 11.5]. Points were slightly shifted 
horizontally by different quantities in order to avoid overlapping. 
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Fig. 21. — M* vs . {U — V) rest frame for our sample galaxies. Solid lines separate the red 
and blue sequences at each redshift (from Bell et al. 2007). 



This figure "f4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "f5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "f6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "fl2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "fl3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "fl5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "fl7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "fl9.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



This figure "f21.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0707.3415vl 



